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TEST PARTICLE METHOD IN KINETIC THEORY 

OF A PLASMA 

Kyoko Matsuda* 

Laboratory for Theoretical Studies 
ABSTRACT 

A general relationship has been found between many body correlation func- 
tions and the conditional probability functions which describe the shield clouds 
surrounding test particles. This is the generalization of Rostoker’s superposi- 
tion principle. The relationship is useful because the problem of kinetic theory 
is reduced to determining the conditional probability functions, wliich involves 
essentially only the Vlasov equation. This method has been applied to obtain a 
second order correction to the Lenar d- Bale scu collison integral. 


*NRC — NASA Resident Research Associate. 
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TEST PARTICLE METHOD IN KINETIC THEORY 


OP A PLASMA 

INTRODUCTION 

1 

Several years ago, Rostoker found a relationship between the two-particle 
correlation function G(X t X 2 1) and a conditional probability function P(X' | Xt) 
which characterizes a shield cloud of field particles of coordinates X ” (x, v) 
surrounding a test particle of coordinates X' = (x' , v' ) . It is 

G(X 1 X 2 1) = f(X l t)P(X 1 | X 2 t) + f(X 2 t)P(X 2 | Xj t) 

+ n J dX' f(X' t)P(X' | Xj f)P(X' | X, t) . (1) 

This relation has been established without solving explicitly for anything and 
has none of the usual restrictions such as spatial homogeneity, adiabatic time 
behavior etc. usually necessary for obtaining explicit solutions. It is useful be- 
cause the problem of ldnetic theory is reduced to determining P, which involves 
only the Vlasov equation. 

o 

A similar relationship has been obtained by Rostoker and the author to 
higher order in the plasma parameter using the procedure of ensemble averages 
of the singular distribution function for a system of initially uncorrelated 
particles. The first two relations are, to second order in the plasma 
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parameter 


G(X l X 2 t) « n J dX/F(X/)P(X/ | X l t)P(X; | X 2 t) 

+ j 7 j dX.'dXj'FCX.^FCX^QCX/Xj' | X,t)Q(X 1 ' X,' I X,t) (2) 

T( 123) = n J dl' FCl'^Cr | 1) P(l' I 2)P(1' I 3) 

*Z>J dl ; d2'F(l') F(2')P(1' I 1)^(2' | 2)3(1' 2' | 3) (3) 

p 

where 

G( 12) - A(l,2) f(l) + G( 12) 

T( 123) = A( 1, 2) A( 2, 3) f ( 1) + A(l,2) G(23) + A(2,3) G(31) 

+ A( 3, l") G( 12) + T( 123) 

A( 1 , 2 ) - ^-8 (Xf * 2 ) 8 (y 1 -v 2 ) 

and G(Xj X 2 t)is written as G( 12). G and T are the usual two-particle and 
three-particle correlation functions, and F is the one-particle function at 
the initial time. £( 1 ' | 1) is the probability of finding a particle at ~ (\ t , v ) 


2 


at time t given a particle at X/ “ (Xj'.Vj* ) at the initial time. jp( 1 ' | l) 

+ P( 2' | 1) + Q(l' 2' I 1) is the probability of finding a particle at (X { t), 
given a pair of particles at (X/ X 2 ' t " 0) . The equations for ?* and Q are ob- 
tained by carrying out averages on the Klimontovitch equation, 

From the definition and also Equations (2) and (3), it is expected that 
P( 1' | 1) must have singularities. For instance, consider a uniform field-tree 
system. Solving the equation of P, we obtain 

P(l' | 1) s ^(Vj - Vj' ) S(x 1 -Xj' - Vj t) + nonsingular terms . (4) 

The first term is interpreted as the probability of finding the test particle while 
the second represents a shielding particle. However this singular solution may 
be criticized because the expansion itself breaks down in the vicinity of test 
particles. In fact, using this solution Equations (2) and (3) do not give the proper 
singularities of G and f . A further criticism of the singular function of Equation 
(4) is given in Appendix. 

Since the main object of Reference 2 is to make the relationship clear be- 
tween kinetic theory and weak turbulence theory, the inconsistency of the solu- 
tion P does not affect the discussion. A correction should be made to the dis- 
cussion concerned with discrete particle effects. Of course, by the termination 
of the hierarchy, the approach in Reference 1 also breaks down in the vicinity 
of the test particle. In both approaches, we are not free from the problem of 
cut-off at short distances. However this is not our present problem. The ap- 
proach in Reference 1 is selfconsistent. The difficulty in Reference 2 is that due 
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to the inconsistency, the restriction to a system of initially uncorrelated particles 
is not able to be removed. Besides the physical pictures of i* i 1), 

$(1' 2' i l), • • * are less clear than that of P( 1' [ 1). 

In this paper, the direct generalization of Reference 1 to higher order in the 
plasma parameter is presented. The general relation as well as Equation (1) has 
none of the usual restrictions, which we have stated below Equation (1), Ros- 
toker’s superposition principle, that is, Equation (1) has been often mi sunder ~ 
stood at this point . Through this relation, the Bogoliubov -Born-Green- 
Kirkwood-' Yvon hierarchy is reduced to the hierarchy of the test particle func- 
tions which involves essentially only the Vlasov equation. 

As an application of this method, a second order correction to the Lenard- 
Balescu collision term has been calculated. The correction is significant for a 
system with a one dimensional velocity distribution such as a plasma with an 
extremely strong magnetic field, for which the Lenard-Balescu term vanishes. 

KINETIC THEORY 

We consider a gas of charged particles interacting only through Coulomb 
forces. The system, may be described by the B-B-G-K-Y hierachy. If we 
introduce an additional charged particle into the system, the field particles form 
a cloud surrounding tfye test particle. To first order in the plasma parameter 
the cloud is characterized by the perturbation P( X 1 / | X x t) of the one-particle 
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distribution function f(X t t) which is govorned by ] 


fd \ ^2 (9f(Xj t) ^ i 

ijt 4 0 , CX 1 't) + 0(X 1 t))p(X l Mx i t) ' ^ T' < 5 > 


where 


0 (X/ t) y‘ ' m CXj 1 t) * 


and 


0(X 1 t) 


V. * 7 


<9 e 


<9 


l ' 0x t m ‘ dVj 


ne^ g££V) £ 
m dY l dx x 



F M (Xt) is the macroscopic force, that is, the external electromagnetic force 
plus the macroscopic Coulomb force 


d 

no 



f(X 2 t) t 
|x -% 2 | 


Infinite mass randomly distributed ions are assumed since the generalization 
including ion distributions is trivial. In Reference 1, the relation expressed by 
Equation (1) has been proved between the solution of Equation (5) and the two- 
particle correlation function. 

In second order theory, we are concerned with a new quantity W(Xj' | X 1 X 2 t) 
which is the perturbation of the two-particle correlation function G(X t X 2 1) 
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duo the the tost particle. To save space, wo introduce singular functions 


P(i‘ 1 1) " P(i' ! i) + &a',i) , 


*( 1' I 12) " W(l' I 12) f A( 1, 2) Pf 1 ' | 1) . 


P(l' I 1) and W(l' I 12) obey the equations 

{&+ O'(l') + 0(1)} <*(!' | 1) 


nc 2 dP(l' 1 1) 
*~ m dY x 


no 2 dPCV \ 1) 
m ” 



P(l' | 2) 


P(l' I 2) 



I 12 > 

■Tv7 W'<1' I 12) 
1 


( 6 ) 
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+ O'(l') +0(1) + 0(2)} W(1‘ | 12) 


~ lie2 
m 


+ 



+ 



2 



J 


e?( 


era 


ere 


Aiin . r 

5V/ J 


d3 


d 

1 

3*7 

| X 1 ~ X 3| 

a 

1 

**i' 

l x ,' ~ X 3 I 

a 

1 

dX 2 

| X 2 ~ X 3 | ' 

d 

1 


G(23) 


G( 13) 


5x/ \x{ -x| 


G( 13) 


nef 5G( 12) f a 1 

m dy % J 0,5 5x t |x t -* 3 | 

ssl i&i2> f a i_ 

n> dv 2 J 0,5 |* 2 -* 3 


P(l' I 3) 


P(l' | 3) . 


(7) 


In the right hand side of Equation (6), the first two terms are effects of Coulomb 
forces due to the dressed test particle; the third and fourth term are perturba- 
tions of the collision integral due to the test particle. The function W' is not 


obvious at this stage. 


The solution of Equation (7) may be given by making use of Equation (X). 
Using the singular functions Equation (1) is written as 


6(12) - n J d2' f(2')P(2' | 1)P(2' | 2) . 


( 8 ) 


The test particle at X t ' produces the perturbation 


W( 1 ' I 12) 


n j d2' 


P(l' I 2') P(2' | 1) P(2' I 2) 




-| d2 ' 


f(2') Q( 1' 2' | f) P(2' | 2) 


+ 


n Jd2' 


f(2') P(2' | 1) Q( 1 ' 2' I 2) . 


(9) 


Q(1'2 I 1) describes the cloud surrounding a pair of particles, in other 
words, the perturbation of f(X x t)due to a pair of test particles at X/ and X 2 ' , 
and obeys the equation 


+ o'(l') + 0'(2' ) + 0 ( 1 )} <X 1 ' 2' | 1) 


ne 2 dP( 
___ ^ 


ill i>. f 

5v i J 


^ <9x t x t - x ? | P(2' | 2) 


ne 2 <9P( 


m 


V 1 1) f 

d \ / J 


d2 


dx/ x; »x 2 


P(2' | 2) 
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M r!q I >»’ l2> 


, ^ rfeUi, . f d2 ^ f 

m dv 2 ' J dx 2 ' |.\ 2 ' “X 2 | ' ' 


This equation is much easier to solve than Equation (7), and the quantity Q has 
much simpler meaning than the quantity W does. 

It is straightforward although tedious to show that relationships exist be- 
tween the first two correlation functions of the hierarchy and the functions P and 
Q, provided that the quantity W' is given by 




W'(l' I 12) = W(l' | 12) - n d2' f(2') P(2' | 1) Q( 1 ' 2' 2) ! (11) 


that is, W' ( 1' | 12) does not contain the term of W( 1' | 12) in which the test 
particle interacts only with the shielding particle. The relationships are 


G(12) 


= njdl' 


f(l') P(l' | 1) P(l' | 2) 


"■J 


dl'd2' P(l' | 2') P(2' | 1') P(T I 1) P(2' | 2) 
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n 2 f 

~2 dl' d2' f(l') f ( 2 ' ) Q( 1 ' 2' | 1)Q(1' 2' I 2) 


+ 2"] n 2 f dl' d2' f(l') P(l' | 2') P(2' | 1) Q(i' 2' | 2) (12) 

P *' 


T( 123) = n I" dl' f(l') P(l' | 1) PC 1 ' I 2) P(l' | 3) 

+ Y"’ n 2 | dl' d2' f(l') P(l' | 2') P(l' | 1) P(2' | 2) P(2' | 3) 

P ^ 

+ ZZ n2 l dl/ d2 ‘ f(1 f < 2 '> I V P( 2 ' I 2) QC !' 2' I 3) (13) 


where 


L 

P 

means the sum over all cyclic permutations of particles involved. Each term 
in Equations (12) and (13) is shown schematically in Figures I and II in the same 
order. 

The rule governing the general relationship between any higher order cor- 
relation function and the test particle function P( 1' | 1), Q( 1 ' 2' | 1), 

R(l' 2' 3' I 1) etc. can be obtained from these diagrams. A dot indicates a test 


10 



particle and a circle indicates a shielding cloud. Each dot counts as order 1A, 
each circle surrounding one center counts as order c , each circle surrounding 
two centers counts as order e 2 , etc. , where a is the plasma parameter. In each 
diagram, every center is surrounded by the same number of circles and the 
total order of a diagram is given by the product of orders for circles and dots 
in the diagram. It should be noted that the interaction between additional par- 
ticles which will be averaged always comes in as a nonsingular function. 

If we rewrite the relationship in terms of nonsingular functions, for the two- 
particle correlation function 


G( 12) = f(l) P(1 | 2) + f(2) P<2 | 1) + n J dl' f(l') P(l' | 1) P(l' I 2) 


ZM 


+ P(1 j 2) P(2 j 1) + ^ n | dl' P(1 | 1') P(l' | 1) P(l' |2) 

p 


+ 


n 2 f dl' d2'P(l' | 2') P(2’ | 1') P(l' I 1) P(2' | 2) 


ZM-' 

l) " 


f(l') PCX' | 1) Q(l' 1 | 2) 




d2' f(l') f( 2' ) Q( 1' 2' 1) Q(l' 2 ' | 2) 


+ 


ZZ-’J" 1, 


d2' f(l') P(l' I 2') P(2' l 1) Q(T 2 ' | 2) . 


(14) 
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This is schematically expressed in Figure in and the three particle correlation 
is shown in Figure IV. The rule governing these diagrams is the same as before. 
The differences are as follows. An open dot indicates one of particles involved 
which itself acts as a test particle and has order e°. There is no restriction 
on the number of circles surrounding this kind of dot. 

The relationship always exists provided it is satisfied at the initial time. 

Now kinetic theory is reduced to the determination of P(l' | 1), Q( 1' 2' I 1) , 

1^1' 2' 3' | 1) , etc. The equations for these quantities are essentially the Vlasov 
equation. 

Rostoker has shown that the first order superposition principle gives the 

known two-particle correlation function for an equilibrium system. In the 

second order theory, we also obtain the equilibrium two- and three-particle 

4 

correlation functions which have already been obtained by other methods . The 
entire procedure is considerably more complicated than that in the first order 
theory and the details will be omitted in this paper. However it may be worthy 
to notice that we obtain some integrals which are hardly calculated by a straight- 
forward way. They are, for instance 



F(u) 

e(k, -iku) | 2 


1 

e(k, 0) ’ 


1 

^ tl 0) 





1 

“i r) e(k 3 , 
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+ e(k 2 


feL dr - im { 


dr r Im ]e(|< 3l -ir) c(k 




* e(k v O) 


£ dT r Im {e(k lf “iT) e(k 2 , 


K* K* K* 2 H 

+ H* + 


_ rjj 


kj 2 k 2 2 k 2 2 k 3 2 k 3 2 k x 2 k l 2 k 2 2 k 3 2 

cCkj.O) e(k 2 ,0) fe(k 3 , 0) 


for finite positive k lt k 2 , k 3 , where F(u) is the one dimensional maxwell function 


F(u) 


u 


& v h, exp r • 


k = w p /v th , and 


e(k, ”ir) = 1 + 


K * 

k 2 


f du u F(u) 

J u - | § 


They are obtained by changing order of integrations in 


J dv i dv 2 G k < V l> V 2> 


or 


J dv i dv 2 dv 3 T k,k 2 ( V l. V V 3> 
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where G k and T k jk ^ are the Fourier transform of correlation functions and ex- 
pressed in terms of the test particle functions, 

COLLISION INTEGRAL 

In this section, we calculate a second order correction to the Lenard- 
Balescu collision term. To this end adiabatic time behaviour and homogeneity 
of the system are assumed. 

The basic technique used here is the Fourier- Laplace transformation. 

From Equation (5), the asymptotic solution P k ( 1' | 1) is easily obtained to first 
order in the plasma parameter. Incidentally making use of Equation (1), this 
solution gives the Lenard-Balescu term. When we know the right hand side of 
Equation (10) and note that 


/ 


dv 2 P k (1' 


2 ) 


1 1 

n e(k, -ik • v i / ) 


(15) 


where 


e(k, -ik • Vj' ) 



dv ik • 0y f 
ik • (v - v ' ) + b 


and damping terms are neglected, we may obtain ^ k (1' 2' | 1). The right 
hand side of Equation (6) is now known to obtain the second order solution for 
P( V | 1). We note that we may use only time independent terms of the asymptotic 
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t 

solutions for the right hand sides of Equation (6) and (10), since all others are 
damped after the integrations over velocities of shielding particles. 

We thus obtain 



(12) dv 2 



f(2)P.^ ) (2[ 1) 
e(k, -ik • v 2 ) 


n] dv,' f(l')P-kV I 1) Jdv 2 P k (2) (l' I 2 ) 


♦ f-*' 


J (2")‘ 


! 1) J 


Pk'k' (!' |2)P k S U C2|l l ) 
2 e(k, ~ik • v 2 ) 


CD 


dv. 


+ n J dv/ f(l')J dv 2 — 


e(k, -ik * v ) 


+ 


n 2 J d»,' dv 2 ' f(l') P_ k i k < 


(1' I 2')P. < k 1> (2' I 1) J dv 2 Q li+k - 2' \2) 


T j dVj'dVj' f(l') f(2') Q-k-k ' , h ' Cl' 2' ! 1) f dv 2 Qk +k < ,_ k '(l' 2' | 2) 


(16) 


where P k (1 ) ( 1 ' | 1 ) , P k (2 ) ( 1 ' I 1) and ^ k ^ ( 1 ' 2 ' | 1) are the time inde- 
pendent parts of the corresponding asymptotic solutions 


P k C1> (l' I 1) 


OJ T 


ik d f 

k 2 ‘<? v 


n e(k, -ik • v') ik * (Vj -v 1 ' ) + S 


(17) 
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Qk,,k s Cl' 2' I 1) = ^Uk, +k 2 , v,: k, •r,'+k.- vj) 


ik 


x A ^ i 3 / j \ 

7(k v (2 ' h) 


ik 




c(k 


4 1 p(o | u 

e(k 2 , -ik 2 ■ v 2 > k,* av i k i +h 2 (1 l 1 ) 


ik 


— pj” (2'|l) 


c(k,, -ikj k,2 r)V 2 ' k l" i 2 


(18) 


d < 2 >/i> ! n - JLi/i i , f elk' ik' 
P k (1 |1) - n L(k,r 1 ik-v^J — — 


e(k', -ik' • Vj') ^ v i 


- p k - k ' (i' li)+-r; 1 


e(k', -ik' • Vj) dr l 


; iw (i' 1 1) 


+ ^?k (i' 1 1) 


(s(k\ -ik' • Vj ) ' X )| 


+ -2- / 

+ dv ^ 


[ r P k (1) (l' |2)P k S 1) (2| 1) r 

L J s (k , -ik' • v 2 ) J 


dv. 


f (2)Q k ,-k' a' 2 1 1) 

e(k', ~ik' *v ) 


n 2 Jdv 2 ' f^2' 


f l2') Pk-k' <2' I 1) 


J *2 Qi,k'- 


k (1' 2' 2) 
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*V*' ( 1 ' I 2 ) P|( ( 2 i 1 ) f f( 2 ) QU'-k' O' 2 I 1 ) 

— — \ n dv 0 — — rT TT. — : — — 


C(k', -ik' *v ) 


2 * (k'» -ik' * v ) 


The operator L is defined as 


L(k, Vjj k • v/) ^(Vj) 


2 ik df 

P k 2 ^ Vl C dv tfC v) 

ik - (Vj *“V j' ) + b K 1 f.(k, -ik * v^) J ik * (v -v/) * 


The kinetic equation is written as 


ii 

at * 




where J < 1 > is the Lenard-Balescu collision integral and 




dk k 

(2jt) z k 2 


dv 0 G/ ; ( 12) 


The imaginary part of Equation (16) can not be expressed in a simpler way un- 
less the one-particle function t‘ is specified. In a system of one dimensional 
velocity distribution where 


Im ‘ J dV 2 °k l> < V 1 V 2> 
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vanishes, we may observe that this second order correction does not vanish in 
general. For this conclusion, it is enough to chock that the term with a particu- 
lar power of functions F(u { ) , 0F/dv \ x , 0 2 FA?u t 2 , • • * does not vanish for 

a given k where F(Uj ) is the one dimensional velocity distribution function, 

For instance, the terms containing <9 4 F/flu 4 , which is the highest derivative 
involved here, result from the first, second and last term of Equation (10), and 
have nonvanishing coefficients. 

DISCUSSION 

The superposition principle reveals the structure of many body correlation 
functions. 

An interesting application of the higher order theory is a systematic study 

c /> 

of mode coupling for an unstable plasma. Recently Price and Harris have 
shown that for the time behavior of the distribution function f , many other effects 
are involved than those considered by quasi-linear theory. Price's work is based 
on Dupree's method and the long time evolution is treated as a succession of 
initial value problems. The limitation of his work is its restriction to a homo- 
geneous plasma and the neglect of three particle effects. The first order test 
particle method gives the same result as Price's for a homogeneous plasma. It 
is expected that a systematic study of three particle effects yields a much different 
result than that given by mode coupling in quasi-linear theory. 
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Since the higher order equations are still complicated, some simplification 
is desirable to apply them to actual problems. To this object, the characteristic 
of each term will be studied in a future paper. 
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APPENDIX 

We show another example in which the singular distribution of Equation (4) 
fails. Using a singular distribution function 


f(X.t) = ±22 S( x ~ x i (t)) 

i 

we have the well known relations 


j d {X i }D N ({X i },t)?(x,{X 1 >) = f(X,t) (A-l) 

J d { X i} D t<({ X i}' t ) f (Xi.{Xi})f(x 2 ,{X i }) = A (1,2) f(X,,t) + f 2 (x t x 2 , t) (A-2) 

etc. , D n is a solution of the Liouville equation* By Liouville’s theorem the left 
hand sides are equivalent to 


J =°) f ( X ’{ X i({ X io}.t)}) 


(A-l') 


j d {*»X({ X i»}'°) f ( X .'{*i ({ X io}- ‘)}) f( X 2 >{ X i ({ X io} ' fc )}) < A - 2 ' 


etc., where X i0 is the initial position of the i tk particle in the phase space. For 
convenience, consider a uniform field-free system. Equation (4) is equivalent 
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to setting 


^,{x i ({X i o},t)}) r n^] 8 ( v “ v io) v(x-x i0 ~ v i0 t) + nonsingular tern . (A-3) 

7 

Although this function is often used in the literature , this function reduces 


(A-2’) etc. to f(X»t s 0) + • • % A (1,2) f (X 4 , t “ 0) + nonsingular 
terms, etc. 

Thus we may claim that the singularity of t or P can not be expressed in 
any simple way with functions of initial coordinates {x. 0 } and time. 
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